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A set of inclusive cross sections at zero rapidity, da/dy\y=o, is presented for pp interactions 
at center of mass energy ^/s = 1.8 TeV. Six elementary particle cross sections are corrected for 
secondary contributions from decays of higher mass resonances in order to produce a set of primary 
cross sections. The primary cross sections per spin state are well described by da^ / dy\y=o ~ 0.721 • 
(tvXI) ■ e-™/^, where m is the particle rest mass, T — hc/r^, and = 0.97 fm. The deuterium 
production cross section is also described if r^ is replaced by = r^A^^^ . The same exponential in 
m and T describes primary charm fractions in e'*'e~ collisions at least up to the J/'tp mass. There 
is no significant evidence for strangeness or charm suppression if only primary production of light 
hadrons is considered. There is evidence that the primary cross section for each particle may have 
the same value for pp and pp collisions and that it may have nearly constant values between -^s — 63 
GeV and y/s — 1800 GeV. Fits to the final state transverse momenta of the particles using a gas 
model favor a temperature T=132 MeV, a chemical potential /j, = 129 MeV, and a transverse flow 
of the gas with /?/ = 0.27. 

PACS numbers: 13.85.Ni,12.40.Ee,13.60.Hb,25.75.Ld,25.75.Nq,25.75.Dw,31.15.Bs,24.10.Ps 



I. INTRODUCTION 

The experiment E735 at Fermilab studied pp colli- 
sions in the Tevatron Collider at center of mass energy 
y/s = 1.8 TeV. The production of particles and their 
pt spectra at rapidity y=0 were recorded along with as- 
sociated multiplicity and any other accessible informa- 
tion. These results have been published in several arti- 
cles P, S i, i, S, i, 0, Hi- A few of the production rates 
were reported as absolute inclusive cross sections, while 
others were reported as rates relative to some simultane- 
ously observed production processes. The purpose of this 
article is to assemble the various rate data and to present 
them all as absolute cross sections. Such a presentation 
reveals patterns which are relevant to specific produc- 
tion models such as the statistical model, string model, 
parton collision models, or the quark gluon plasma. 



II. DATA ANALYSIS 

During the data taking period the accelerator opera- 
tors made periodic measurements of accelerator param- 
eters with flying wires @ and wallcurrent monitors [To| 



to enable a calculation of our luminosity as a function 
of time. A careful analysis of the luminosity data was 
subsequently tabulated by Gelfand, Grosso-Pilcher, and 
White [ll|, [l2| and made available to experimentalists. 

The luminosity data can be used to obtain an absolute 
cross section for the production of any one particle species 
from the relation 



iV = cr 



C{t) dt, 



(1) 
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where £,{t) is the measured luminosity as a function of 
time, cr is some cross section of interest, and N is the 
number of events which are associated with that cross 
section during the time interval T. We have previously 
used the luminosity tabulations together with Eq.([T]) 
to calculate and publish invariantproduction cross sec- 
tions for photons Q, phi mesons deuterium, and tri- 
tium 0. Luminosity issues are discussed in more detail 
later in this article. 

We have repeated our previous analysis of tt, K, and p 
inclusive production with some refinements that were in- 
tended to reduce background, minimize tracking losses, 
and improve resolution in the time-of-flight mass spec- 
tra p^ . The new n/p ratio remained within 1.3% of 
that reported in previous publications, but the relative 
charged kaon flux increased ~1.4 standard deviations 
with the new analysis. Raw data were corrected for de- 
tector and trigger efficiencies, trigger pre-scaling, trigger 
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dead time, and particle loss in 0.09 absorption lengths of 
the spectrometer arm before the last time-of-flight wall 
at 4 meters. Absolute cross sections d(j/dy\y^Q were cal- 
culated for these three positively charged particles. 



III. MEASURED CROSS SECTIONS 

TableUgives a list of measured single particle cross sec- 
tions as collected from various sources. Column 1 lists 
the particles studied. Antiparticle cross sections were the 
same as particle cross sections within errors. However, 
the cross section for antitritium was not measurable. Er- 
rors listed here are the result of convolution of statistical 
and systematic errors. Information on the associated par- 
ticle multiplicity and the invariant pt spectra du/dydpl 
for all these particles can be found in references cited in 
the table. 

The cross section listed for inclusive single photon pro- 
duction is twice the inclusive cross section of a charged 
pion within errors. If one assumes production of a neutral 
pion is equal to production of a charged pion, this value 
is consistent with the assertion in reference [H] that single 
photon production originates with the two-photon decay 
of the neutral pion and possibly some small admixture of 
eta decay. 

In calculating invariant cross sections for the A" and 
S^, we make use of particle ratios in reference Q such 
as [{hP + hP)/{p + p)\ "in the spectrometer arm". Al- 
though it is not explicitly stated in reference [1], these 
actually are rapidity based ratios such as [dN{hP + 
K°)/dy\y=o\/[dN{p + p)/dy\y=Q\, where dN/dy is the 
number of particles per unit rapidity. 

TABLE I: E735 values of da/dy at rapidity j/ = for various 



particles produced in pp collisions at ^/s 


= 1.8 TeV. 


Particle 


da/dy\y^o 


Reference 


7 


114.0 ± 7.8 mb 


[5] 


■K+ 


56.0 ± 6.3 mb 


this paper 


K+ 


7.37 ± 0.79 mb 


this paper 


P 


4.37 ± 0.62 mb 


this paper 


<t> 


0.763 ± 0.202 mb 


[6] 


A° 


1.77 ± 0.41 mb 


[8] 




0.282 ± 0.082 mb 


m 


d+ 


2.02 ± 0.49 nh 


m 


t+ 


0.83 ± 0.34 /xb 


m 



In the literature one frequently encounters experimen- 
tal particle production ratios that appear to be defined 
by the angular aperture of a detector element. For such 
cases the ratios quoted are more likely to be formed with 
dN/dr] rather than with dN/dy, where drj is an inter- 
val of pseudorapidity which depends only on production 
angle and not particle mass. Other discussions of par- 
ticle production ratios make use of the scaling variable 
X = Ipij^fs to compute dN/dx, where pi is the longi- 
tudinal component of momentum in the center of mass 



system. 

In order to compare our data with such measurements 
one can use the approximate conversion factors in Table 
Ullto convert our values of dN/dy to dN/drj or to dN/dx. 
In this table we present conversion factors Cr/ and that 
can be used to multiply our values of dN/dy for an ap- 
proximate conversion to i] or x variables at y=0. 

{dN/dy)y^o = {l/c,j){dN/dr])r^^o = il/c^){dN/dx)^=o 

(2) 

where 

= Pt/ \/p}+m? and c^^ = {^/s/2)/yJp1 -f m?. 

The values of and Cx are approximate to the extent 
that we have assumed pi = and used observed average 
values of pt for each particle mass to estimate in Table [Tl] 
the conversion of dy to equivalent drj and dx intervals in 
our experiment. 

TABLE II: Columns three and four list factors by which one 
can multiply dN/dy in this experiment to obtain approximate 
values of dN/dr] and dN/dx for the particle masses listed in 
column one. The average values of particle pt listed in column 
two were used to evaluate and c^. 



mass (GeV/c^) <pt> (GeV/c) c,, Cx 

7(0) 0.192 ±0.007 LOOO 4688 

7r(0.139) 0.433 ±0.020 0.952 1978 

is:(0.494) 0.590 ±0.020 0.768 1170 

p(0.938) 0.716 ±0.030 0.607 763 

(/)(1.019) 0.94 ±0.26 0.678 649 

A(1.116) 0.750 ±0.025 0.558 670 

H(1.321) 0.90i;',-22 0.56 3 5 63 



IV. CROSS SECTION CHARACTERIZATION 

Figure [T] is a plot of production cross sections versus 
particle mass from Table HI One observes what appears 
to be a rough exponential decrease in cross section as the 
particle mass increases. This decrease is not just some 
overall energy or phase space limitation. At most about 
200 charged particles are observed to be created in a few 
of the highest multiplicity collisions when the available 
energy would allow creation of around 6,000 pions, each 
with a kinetic energy of a pion rest mass. Most of the non 
diffractive particle production is far from encountering a 
simple phase space limitation. 

Some features of Fig. [T] are noteworthy. The photon 
cross section is not an independent data point since it 
is mostly attributable to ir^ production Likewise 
electron/positron production is negligible. One can ob- 
serve electrons in the small momentum window between 
50 MeV/c, below which the field of the analysis mag- 
net rejects all charged particles, and 200 MeV/c, above 
which pions begin having enough energy to reach the fi- 
nal time-of-flight counters. The electron mass peak is 
sharply defined, but it disappears completely above 250 
MeV/c underneath the low mass tail of the pion mass 
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FIG. 1: Cross sections da/dy\y=o from Table |T] plotted as a 
function of particle mass. The dashed line is an arbitrary 
exponential falloff inserted to guide the eye. 



peak. The observed part of the electron flux is consistent 
with photon conversions in the 2 mm thick aluminum 
wall of the beam pipe. Since electrons and photons do 
not participate in the exponentially falling behavior of 
the cross sections, we will omit them from most of the 
remaining discussions. 

Although the two nuclear cross sections (d and t) are 
small, as the trend for elementary particles suggests, their 
behavior may involve unique nuclear effects. We will de- 
fer treatment of the nuclear cross sections in a quantita- 
tive way for a later section. 



pression, all cross sections would be described by 



da 
dy 



(3) 



where m is the particle rest mass. We then explored 
variations of strangeness suppression for the following 3 
cases: (a) A single divisor A raises all strange particle 
cross sections equally, (b) A single divisor A is applied to 
measured cross sections for and A" because they con- 
tain only one strange quark, but a divisor of A^ is applied 
to 4> and S~ cross sections because each of these parti- 
cles contains 2 quarks with strangeness quantum num- 
bers, (c) A divisor of A is apphed to , A°, and 4> cross 
sections because only one quark pair must be produced 
to form these particles, but a divisor of A^ is applied to 
the S~ cross section because two quark pairs must be 
produced to form it. 

The data are fitted for values of A, B, and A that min- 
imize x^. The values of per degree of freedom for 
each of the cases listed above are: (a) Xdo/=4.43, (b) 
Xrfo/=l-33, (c) Xrfo/=4.37. The best of these fits is case 
(b), and the dashed line in Fig. [2] shows graphically the 
extent to which the adjusted strange particle cross sec- 
tions fit for A = 0.437. 



P0 A 



^SIGMA=80.0*EXP(-3.02«M) 



V. STRANGENESS SUPPRESSION MODELS 

Most strange particle production in Fig. [1] appears sys- 
tematically suppressed relative to non-strange produc- 
tion even after allowing for some exponential decrease 
of cross sections with mass. Efforts to characterize this 
strangeness suppression have a long experimental and 
theoretical history. A strangeness suppression factor 
A = 0.5 was adopted in the Feynman-Field jet fragmen- 
tation model to address this observation in experimental 
data Other authors have presented evidence for a 

value closer to A = 0.3 

There is probably nothing fundamental that one can 
learn from the measured cross sections in Table |T] without 
correcting them for secondary contributions that come 
from higher mass decays. Nevertheless, we have exam- 
ined the 6 elementary hadron cross sections to see if our 
data might be made to appear more uniform with re- 
spect to mass dependence by using a single strangeness 
suppression factor A. To some extent this effort will serve 
to compare our raw data to that of previous experiments. 
We have explored 3 variations of strangeness suppression. 

We first assumed that if there were no strangeness sup- 
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FIG. 2: Six inclusive hadron cross sections of Table |l] are 
plotted after multiplying and hP cross sections by a factor 
A~^ and multiplying and H~ cross sections by a factor . 
The best fit to an exponential in mass is the dashed curve for 
a strangenss suppression value of A=0.437 with a xlo/ — 1-33. 



If we wish to attach some physical significance to the 
fit such as interpreting the constant B as inverse temper- 
ature, then we should use Eq.(|3]) below instead of Eq.® 



da 
dy 



y=Q 



(4) 



where (2J+1) is the spin degeneracy of a particle or 
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resonance. To see exponential behavior one must plot 
[1/(2 J + l)]da/dy\y^Q. These values of cross section per 
spin state are given in Table Hill 



TABLE III: Values of da/dy at j/ = for hadronic particles 
produced in this experiment by pp collisions at ^/s = 1.8 TeV. 
Spin degeneracy is removed in the last column. 



Particle Mass (GeV) da/dy\y^o 


da/dy\y=o/{2J + 1) 


n+ 


0.139 


56.0 ± 6.3 mb 


56.0 ± 6.3 mb 


K+ 


0.494 


7.37 ± 0.79 mb 


7.37 ± 0.79 mb 


P 


0.938 


4.37 ± 0.62 mb 


2.19 ± 0.31 mb 


<P 


1.019 


0.763 ± 0.202 mb 


0.254 ± 0.067 mb 


A° 


1.116 


1.77 ± 0.41 mb 


0.885 ± 0.203 mb 




1.321 


0.282 ± 0.082 mb 


0.141 ± 0.041 mb 


d+ 


1.871 


2.02 ± 0.49 ^b 


0.673 ± 0.163 ^b 


t+ 


2.810 


0.83 ± 0.34 fih 


0.415 ± 0.170 ^b 



The three cases for a strangeness suppression A that 
were discussed above were also fit to the spin-scaled cross 
sections of Table IIIIl The values of per degree of 
freedom for these three cases were: (a) Xdof ~ 10.7, (b) 
Xdof = 3-08, (c) Xdof — 12.06. The best fit is once again 
for case (b) , but the value of Xdof is even less satisfactory 
than the fit to raw data. For this case A = 0.467, and 
interpretation of the mass dependence implies T — 248 ± 
11 MeV. 



VI. THERMAL DESCRIPTION INCLUDING 
PARTICLE DECAYS 

The assumption of a strangeness suppression factor A 
improves the exponential appearance of the data, but it 
does not return a statistically satisfactory description of 
the measured data. Also it is not clear if A results from 
particle production or particle decay or some combina- 
tion of the two processes. 

It is most certain that particle and resonance decays fa- 
vor pion and proton final states over strange particle final 
states. Every particle and resonance ultimately results in 
pions and nucleons, whereas only flavored particles and 
flavored resonances can decay into strange particles. 

We have investigated a thermal model which has no 
strangeness bias in the production of particles and res- 
onances. It does, however, use the actual strangeness 
bias found in decays as tabulated by the PDG ^IS^ for 99 
of the lowest mass particles and resonances. The model 
assumes that a particle of given mass is produced in pp 
collisions with equal probability in (2 J -I- 1) spin states, 
{21 + 1) isospin states, and particle/antiparticle states. 
Each spin state for a specified charge is assigned a pri- 
mary production cross section per spin state according 
to its rest mass m as in Eq.([5|). 

^ = Ae-"/^ (5) 
dy y=o 



Observed particles are considered to result from either 
primary or secondary production. Primary production 
cross sections , da^ /dy\y=Q, must be estimated using the 
primary particle mass in Eq. ([5]) since experimental values 
for these cross sections are unknown. Secondary contri- 
butions are determined by the primary production rates 
for some higher mass resonances as given in Eq.([5|) fol- 
lowed by a variety of branching ratios in their ultimate 
decay to a primary particle of interest. For low mass 
particles these branching ratios are very well known and 
have errors that range from a fraction of a percent to 
several percent [l^] . 

We are interested in the primary production of a parti- 
cle because only that rate might be expected to follow the 
exponential of Eq. ^ . Our goal is to learn what fraction 
X of the measured inclusive cross section for a particle 
is the result of primary production. Measured cross sec- 
tions da/dy can then be related to primary cross sections 
per spin state da^ /dy by Eq.®. 

,„ _ , , daP da 

dy y=o dy y=o 

The value we find for the fraction a; of a particular 
particle is a function of the assumed temperature T in 
Eq.®. For example, when we try to determine the pri- 
mary 7r+ inclusive spectrum, we may find that a massive 
resonance with large J and / values produces many 7r+ 
mesons by its decay chains, but because its resonance 
mass is much greater than T, it makes a relatively small 
secondary contribution to the observed ■k'^ inclusive spec- 
trum. 

Branching ratios, conservation laws, and Clebsch- 
Gordon coefficients were used to determine how much 
each primary state contributed to observed 7r+, , p, 
(/), A°, and S~ inclusive spectra in the E735 detector. 
For particles with lifetimes of « 10~^° seconds it was 
necessary to use a Monte Carlo program to evaluate av- 
erage corrections for loss of particles due to detector aper- 
ture and analysis cuts. Typically ~65% of these protons 
were accepted and ~15% of these pions. Strong decays 
give prompt secondary particles at the beamline and were 
given no relative corrections for the detector losses. 

The tableau in Fig. O shows, as a function of the mass 
of a particle/resonance, the relative integrated contribu- 
tions from various resonances and particles to the 6 in- 
clusive cross sections we are studying. Vertical scales in 
the figure were chosen arbitrarily in order to illustrate 
better the limiting contributions from decays of higher 
masses. The exact shape and limiting values of these 
curves will have a dependence on the choice of T used 
in Eq. ^ . The value of T was chosen to be 204 MeV in 
order to produce the set of curves shown. The ratio of 
the lowest mass step (primary production) to the limit- 
ing integrated value yields x, the fraction that is primary 
production. 



5 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 2.5 
MASS (GeV/c^) MASS (GeV/c^) 



FIG. 3: Six plots are shown for a common assumed value of T. 
They display relative integrated cross section for a particle as 
a function of the mass of contributing resonance or particle 
decays. In general the limiting values of these curves will 
depend on the choice of T. These plots use the value T — 
204 MeV, which minimizes the chi-square fit of primary cross 
sections per spin state to Eq.©. 

A two parameter fit was done for A and T in Eq.([5]) 
using the 6 inclusive cross sections per spin state in col- 
umn 4 of Table Hm In order to do this it was necessary to 
make the functional dependence of the primary fraction, 
X, on T available to the fitting program. The fit is 
excellent with a P^r degree of freedom Xdof ~ 0.317 
for 4 degrees of freedom. 

Table [TVl lists for each measured particle the value of 
X and the primary cross section per spin state as deter- 
mined by the fit. The primary cross sections per spin 
state are plotted in Figure H] along with the fitted func- 
tion (dashed curve). The temperature determined by 
the fit is T = 204 ± 14 MeV. This is significantly higher 
than T = 160 MeV suggested for the bootstrap statisti- 
cal model [l^ and quite close to T = 200 MeV suggested 
for the QGP 

TABLE IV: Values of da /dy at y — for hadronic particles 
produced in pp collisions at ^/s — 1.8 TeV. Column 2 lists 
cross sections with spin degeneracy removed. Column 3 lists 
factors by which column 2 must be multiplied to obtain pri- 
mary cross sections per spin state listed in column 4. 



Particle 


da/dy\y=o/{2J+l) 


X 


dcr''/dy y=o(mb) 


■K+ 


56.0 ± 6.3 mb 


0.433 


24.3 ± 2.7 


K+ 


7.37 ± 0.79 mb 


0.522 


3.85 ± 0.41 


P 


2.19 ± 0.31 mb 


0.223 


0.489 ± 0.069 


4> 


0.254 ± 0.067 mb 


0.999 


0.254 ± 0.067 


A° 


0.885 ± 0.203 mb 


0.176 


0.156 ± 0.036 




0.141 ± 0.041 mb 


0.527 


0.0743 ± 0.0216 



The extraordinarily good fit to a simple exponential in 



mass is somewhat of a surprise since it relies in no way on 
the significant amount of translational energy observed 
experimentally for these particles. There is also no hint of 
a distinction between Bose-Einstein statistics and Fermi- 
Dirac statistics. We do have good evidence that there is 
some Bose-Einstein interference between observed pion 
pairs in this experiment 2T]. 

The best fit values for A (45.3±2.8 mb) and T (204±14 
MeV) are remarkable for their apparent relation to other 
meaningful physical parameters. As we note in Section 
XV , if cross sections fall exponentially, the Uncertainty 
Principle suggests that T = Tic/rh^ where the hadronic 
radius is ru « 0.97±0.07 fm. The value of ^ = 45.3 mb « 
0.721 • (ttA^) should be compared to the fact that the 
TT^p cross sections [Tsl at the J=3/2, 1=3/2 resonance 
saturate for a value given by 

a(7rp) w 0.797-(7rA^)(2J+l). (7) 
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= 45.3«EXP(-4.90«m) 
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FIG. 4: Measured cross sections are scaled down to remove 
secondary contributions from decay of higher mass resonances 
and then divided by spin degeneracy (2 J -I- 1) before plotting. 
The dotted line is a Boltzmann distribution in mass using 
T = 204 Mev obtained from a fit to the primary cross sections 
for the 6 elementary particles. The two nuclear cross sections 
per spin state are plotted to show the extent to which they 
do not appear to participate directly in this simple thermal 
model. 



VII. EXTRAPOLATION TO HIGHER MASSES 

If particle production as described by the above ther- 
mal model is valid for higher particle masses, then the 
thermal model can serve as a lower limit for presently 
unmeasured cross sections or it can be used as a standard 
measure for identifying unusually large cross sections. 

Checking the validity of the thermal model at high 
mass values is difficult for two main reasons. One rea- 
son is that almost all event triggers at the large collider 
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detectors rely on finding particles with high pt . The few 
cross sections which have been measured apply only to 
the highest portion of some incompletely known momen- 
tum spectrum. 

Another reason is that detailed branching ratios like 
those we used to prepare Fig. [3] are unknown or unpub- 
lished, so that one cannot easily calculate a primary cross 
section. Only if one could find a particle like the 4> me- 
son (Fig. 3d) with negligible contributions from nearby 
mass resonances, would it be unnecessary to know de- 
tailed branching ratios for that particle. 

We have attempted to make ordcr-of-magnitude cross 
section estimates per spin state for four massive particles 
using reported data. These are plotted as X's without 
error bars in Fig. [5] for comparison with the extrapolated 
thermal model. Some calculational details of these esti- 
mates are outlined in the Appendix. 
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However, since the number of secondary sources tends to 
increase exponentially with mass p^ . it may never be 
practical to isolate primary cross sections for the heavier 
particles. 

The estimated J/tp cross section is approximately two 
orders of magnitude above the thermal model prediction. 
There are 11 spin states of Xci, Xc2, and ^{2S) which 
have significant decays (20% to 56%) into J/ti^- How- 
ever, it seems unlikely that the primary thermal produc- 
tion of J /if) is a source which is competitive with other 
constituent scattering sources. 

Momentum distributions have been measured for D*^ , 
J/'f/', and B'^ production perpendicular to the colliding 
beams [H [H, HI]. In the range of 6 < pt < 20 GeV/c 
the D*^ and J /ip spectra fall off with a definite inverse 
power law {du/dpt ~ 1/Pt)- This is the same large pt 
dependence observed in pion production at 90° for large 
^/s collisions (1^ . Independence of ^/s and particle mass 
are sufficient reasons to identify this common behavior as 
Feynman-scaling [23|, which is typically associated with 
hard parton scattering [2^. However, the pt dependence 
of the cross sections we have measured in this experi- 
ment is better represented by an exponential function up 
to 1.5 GeV/c for identified protons and up to 3.0 GeV/c 
for A°'s. At large particle masses where the thermal cross 
sections are expected to be very small, hard parton scat- 
tering may be the dominant production mechanism. 



VIII. HADRONS FROM e+e" COLLISIONS 



PARTICLE MASS iGe^/c") 

FIG. 5: The dashed curve is a fit to cross sections observed 
in this experiment. In the figure the curve is extrapolated 
to much higher mass values and is bounded by dotted curves 
which represent statistical errors for the fitted temperature. 

An extrapolation of the fit to our six primary cross 
sections is shown as a dashed line in Fig. [5] bounded by 
dotted lines representing the slope errors for the fit. Al- 
though the higher mass cross section estimates in this 
figure have varying degrees of unreliability, they all lie 
above the expectations of thermal production, and the 
discrepancy tends to grow with increasing mass. 

The Z?*+ cross section appears to be almost an order 
of magnitude above the thermal production expectations. 
There are very few examples of secondary decay contri- 
butions from nearb y re sonance decays into D*^ that are 
listed by the PDG [18|. Admixtures of B mesons are ex- 
pected to decay as D*+ + Anything 22.5% of 
the time [l^, but the cross section alone is about an 
order of magnitude below the D*^ cross section. 

If all secondary sources of D*+ can be accurately ac- 
counted for, it may eventually be possible to calculate a 
primary cross section that agrees with the thermal model. 



Although reliable cross section values for higher mass 
particles are not available for pp collisions, some have 
been measured in e^e~ collisions. It would be interest- 
ing to learn if the thermal model in some way mediates 
hadron production in electron collisions. We have cho- 
sen to investigate this possibility by studying multiplicity 
fractions for charm particle production. 

The multiplicity fraction fh for a particle is the ra- 
tio of the inclusive cross section for producing the par- 
ticle to the total cross section for producing hadrons: 

= (ji-nci I ^hadron ^ Thus wc treat multiphcity fractions 
for charm production in e^e~ collisions in the same way 
as we treat inclusive charm cross sections resulting from 
hadron collisions. For example, the multiplicity fraction 
for production, fh{D~^), contains a primary compo- 
nent as well as secondary contributions from the decay 
of higher mass particles. 

The PDG provides tables ^] of the total hadronic 
fractions fh for a set of charm particles. Although data 
points are tabulated for collisions at y/s — 91 GeV and 
y/s = 10 GeV, we will focus on the data from ^/s = 91 
GeV in an effort to achieve smaller errors. Values of fh 
measured at -y/s = 91 GeV are listed for 6 charm particles 
in Table El and plotted in Fig.[6ta). 

Because the production rates of b and c quarks in e+e" 
scattering are comparable in size at ^/s = 91 GeV (the 
mass) , the contribution of b quarks to the total hadronic 
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fraction of a charm particle can be substantial and must 
be subtracted in order to obtain the fraction contributed 
by c quarks, 



fc — fh — fb- 



(8) 



Fortunately it is possible to obtain for a limited 
number of charm particles from a set of partial widths 
r tabulated by the PDG. [30(1 Using the partial width 
r(6 — + c) for a b-admixture to yield a charm particle c, 
one can calculate ft directly. 



/b = T{b ^ c) X T{Z" bb)/T{Z° 



hadrons) 



/b T{b ^ c) X (0.2163 ± 0.0007), 



(9) 



(10) 



where T{Z^ bb) and T{Z'^ — > hadrons) are the partial 
widths for Z^ decay into bb and hadrons respectively [3l| . 

The first five values of fc listed in Table |V] were ob- 
tained as outlined above. Errors for these values range 
from 5% to 30%. However, the last value, which is for 
J/ip production, presents a special case because 95% of 
J/ip particles produced at ^/s = 91 GeV come from the 
decay of B hadrons [12] ■ The error on the difference in 
Eq.® may be comparable to the difference itself. 

One might consider associating fc with the experimen- 
tally prompt component of J/ip production at y/s = 
91 GeV, but the statistical accuracy on a direct mea- 
surement of prompt J/ip production at this energy is 
poor [l^]. Nevertheless, the prompt value at y/s = 91 
GeV is close to the total value [2^ obtained with much 
greater accuracy at y/s — 10 GeV (where b quarks are 
not produced). We have used this more accurate value 
for fc in Table IVl for J/ip production. 



TABLE V; This table lists multiplicity fractions extracted 
from charm production in e^e~ collisions [2^. Column 1 
identifies the charm particle. Column 2 gives the experi- 
mentally observed fractions for charm particle production at 
^/s = 91 GeV. Column 3 gives the fraction when contributions 
from b decays are removed. Column 4 gives the primary frac- 
tion after contributions from higher mass charm decays are 
also removed. 



Name 

D+ 
Dt 
D*+ 
At 



fc 



0.454 ± 0.030 
0.175 ± 0.016 
0.131 ± 0.021 
0.194 ± 0.006 
0.078 ± 0.017 



0.325 ± 0.031 
0.125 ± 0.016 
0.099 ± 0.022 
0.156 ± 0.007 
0.057 ± 0.018 



0.0690 ± 0.0065 
0.058 ± 0.008 
0.033 ± 0.007 
0.0384 ± 0.0017 
0.0057 ± 0.0018 



0.0052 ± 0.0004 0.0005 ± 0.00005 0.000148 ± 0.000015 



Our objective is to determine if the 6 observed values 
of fh in Table |V] are consistent with the thermal model 
that describes hadron production in pp collisions. To do 
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FIG. 6: (a) Observed hadronic multiplicity fractions at 
y/s = 91 GeV [2^ vs. charm particle masses, (b) Primary 
multiplicity fractions per spin state after correcting for sec- 
ondary decay contributions. Dashed curve is an exponential 
6.53 X 10^ e'™''^, where m is particle mass and T = 0.204 
GeV. 



this we must correct fc for contributions from secondaries 
of other charm decays. We can write fc as a sum of 
contributions as follows. 



/e=(2J, + l)/,^' + ^(2J, + l)/fB 



(11) 



In this expression fc is the observed multiplicity frac- 
tion for charmed particle c. The primary fraction per 
spin state for particle c is f^, and Jc is the spin of par- 
ticle c. Similarly ff is the primary fraction of particle i 
with spin Ji that contributes secondary decay products 
to the multiplicity fraction observed for particle c. Bic is 
the branching fraction of particle i into particle c, and in 
some cases it is the result of a decay chain. 

In the charm sector the branching fractions Bic are 
sometimes poorly determined. We make use of data 
summarized in 2008 by the PDG [l^. For some parti- 
cles branching fractions into hadrons are given as "seen" 
or "dominant" . In order to obtain numerical estimates 
in these cases we make use of isospin symmetry, phase 
space ratios, and Clebsch-Gordon coefficients. Usually 
the most uncertain branching fractions are for higher 
mass particles which have smaller multiplicity fractions 
and thus contribute fewer secondaries. 

To test for consistency we estimate the primary frag- 
mentation ff using the exponential fall off which we 
found for cross sections in pp collisions. 



(12) 
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where T = 0.204 GeV. After inserting this in Eq.dTTl), we 
can solve for /|', the primary fragmentation fraction per 
spin state of charmed particle c. 

fc 



(2J, 



1) 



(2J, + l)S,,e„ 



(13) 



This expression has two useful attributes. (1) Compu- 
tation of depends only on the measurement of one 
multiplicity fraction, /c, and not on other measured frag- 
mentation fractions and their possible experimental er- 
rors. This could be useful in identifying a single faulty 
measurement since there is no interdependency on other 
multiplicity measurements. (2) The primary fraction 
can never be negative, which can happen if one "sub- 
tracts" secondary contributions using other measured 
fractions that have experimental errors and uncertain 
double counting in decay chains. 



One can conclude that the six charm multiplicity frac- 
tions in Fig. [DJb) are consistent with a thermal produc- 
tion model. We have optimistically assumed that er- 
rors in the branching fractions Bic are negligible. Since 
branching ratios for decays of higher mass particles are 
sometimes poorly determined, these errors in some cases 
could be significant. 

We have also attempted this same analysis using e+e^ 
multiplicity fractions tabulated entirely for « 10.5 
GeV 29'|. By working at that energy experimenters 
should have eliminated b-quark fragmentation from the 
data. Our analysis result at this lower value of ^/s yields 
a fairly good description by T = 204 MeV for the non- 
strange charm particles, but there is a "strangeness sup- 
pression" of about a factor of two for and £'*"'" mesons 
which does not appear at y/s = 91 GeV. This could be a 
physical effect that is related to the choice of a lower ^/s 
or it could be the result of some systematic experimental 
error, such as a slightly incorrect Monte Carlo model. 



TABLE VI: Column 2 of this table lists sources of secondary 
particles that were used to correct the observed fractions for 
the particles in column 1. 



Particle Secondary Sources 

D D* , Do (2400) , Di (2420) , Di (2430) 

Di (2460) , Di (2640) , Dsi (2536) , Ds2 (2573) 
Ds D;(2112), D«o(2317), D,i(2460), Dsi(2536) 

D* _Di(2420),I>i (2430), Da (2460) 

Di(2640),Dsi(2536) 
Ac Ec(2455), Ec(2520), Ac(2595), Ac(2625) 

Ec(2800), Ac(2880) 
J/i^ Xci,Xc2,^(25),A(3872) 



The sources used for secondary decay contributions 
are listed explicitly in Table IVII All of these sources 
have been limited to the charm sector since the earlier 
subtraction of fb accounted for the contributions from B 
hadrons. 

Six calculated primary multiplicity fractions per spin 
state, fP, are listed in column 4 of Table Ivl and plotted 
in Fig. ISlJb). These were computed using Eq. (|13p with 
fc data from column 3 of Table [V] The dashed line is a 
fit for normalization using fP — ^e^™/"'^°^, which gives 
A = 653 ± 24 wfth xlof = 2.8. This could be considered 
a reasonably good fit, since the error bars in Fig. ^h) 
simply employ the same percentages as those given for 
the observed multiplicity fractions, and no allowance is 
made for systematic inconsistencies among various ex- 
periments. 

When corrections for all secondary decay contributions 
are completed, the D^/Z?"*" isospin violation evident in 
Fig. int^a) disappears in Fig. [H^b), which is encouraging 
since we did not use directly measured D* fractions to 
bring about the agreement. This isospin symmetry viola- 
tion (difference in multiplicity fractions) for the and 
is understood in terms of energy violation for decay 
of £>* into D+ 33]. 



IX. MODELS FOR PRIMARY CROSS 
SECTIONS 

The scale factor for our cross sections (ttA^) has been 
a familiar feature since the discovery of the A(3, 3) reso- 
nance, and the temperature (T) appears to be almost 
a trivial consequence of hadronic size and the Uncer- 
tainty Principle. Nevertheless, one would like to find a 
more fundamental description of particle production in 
high energy collisions that could give physical meaning 
to these and other observables. 

In what follows we will examine several descriptions. 
The foremost requirement for each of these is that it must 
be compatible with the observed primary cross sections. 
In addition to that, a useful description should guide us 
to some reasonable approximation for the measured aver- 
age transverse momentum <pt> of each studied particle. 
The measured value of <pt> is, of course, not that of the 
particles from primary production, since there is in gen- 
eral no experimental way to distinguish primaries from 
the secondaries of resonance decays. One can hope that 
the random nature of the decay process will have less ef- 
fect on the average pt than it will on the exact shape of 
the Pt distribution. We present some Monte Carlo calcu- 
lations in Section XI that support this hope. Finally, a 
successful description should yield some gross approxima- 
tion to the observed pt spectra after making reasonable 
allowances for distortions by contributions from secon- 
daries of resonance decays. 

As we investigate descriptions, we will recognize the 
possibility that the final state of the particles may be 
different from the initial production state. That is, we 
consider cases where the final temperature may be differ- 
ent from the initial temperature as long as the particle 
ratios remain the same. 
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A. Hawking-Unruh Radiation 

A recent approach to particle production has been to 
consider it as a consequence of Hawking-Unruh radia- 
tion [s^l • In this description the rate R for the production 
of a particle of mass m is 



(14) 



where a is the particle acceleration. This is the result 
of a quantum mechanical tunneling calculation first used 
by Schwinger to describe pair production by an electric 
field 35j. Hawking later used this method to calculate 
the temperature of black holes, where the acceleration at 
the event horizon was considered to serve exactly as a 
temperature (36l |. 

When applied to high energy experiments, the acceler- 
ation process is initiated by the particle collisions. Some 
estimates for the temperature are T = Qgjl-K k, 200 
MeV m, T V3o-/47r [13, T = ^072^ [13]. Qs is 
the "saturation scale" in the gluon saturation regime. 
In the last two expressions a is the string tension be- 
tween quarks in units of MeV"^ . Our fit to the data with 
T = 204 MeV is in agreement with acceptable values of 
string tension. 

It is important to note that the above estimates of ac- 
celeration do not have a mass dependence. A more recent 
version of this approach [1^ introduces quark masses. 
The resulting temperatures then imply strangeness and 
charm suppression. In earlier sections we have shown 
there is no significant evidence to support suppression of 
strangeness or charm production in primary cross sec- 
tions. 



B. A Bootstrap Description 

The bootstrap model was proposed by Hagedorn [13 
after he first noted that the density of mass states known 
at that time increased almost exponentially with mass. 
A single state is defined as a particle or antiparticle with 
particular eigenvalues for mass, spin (J2), and isospin 
(7z). Frautschi [Isl has written a paper that refines 
the bootstrap assumptions and concisely explains such 
a model and its consequences. The principal results are 
a limiting temperature Tq, sometimes called "the boiling 
point of matter" , and a limiting analytic form for the 
density of states, p{m), as a function of the mass m. 



p{m) = cm" e"/^" (a < -5/2). 



(15) 



Hagedorn has suggested a "smoothed" version of this ex- 
pression which could be used at low masses or even at 
m = (2^,iQ], 

p{m) = cTg/' {ml + m^y^/^ e"/^°. (16) 

Hagedorn expected values of c = 5.5^^^, mo = 500 MeV, 
and To = 160 MeV. 



The inverse of this expression is so similar to the mass 
dependence of da / dy\y=Q that we have tried fitting the 
cross section data with 



daP 
dy 



= AT 



-3/2 



(mi 



2n5/4 ^-m/T 



(17) 



where A and T are fittable parameters, and mo ~ 500 
MeV is Hagedorn's choice of a smoothing parameter. 

The results are shown in Fig. [7] with A = 17.2 ± 1.5 
mb/GeV, T = 0.1397 ± 0.0032 GeV, and xL/ = 0.30. 
The parameters A and mo are strongly correlated and 
poorly constrained by the cross section data, which is 
why we accept Hagedorn's determination for mo with- 
out making it a free fitting parameter. We observe that 
To > T, so that the cross section appears to decrease 
with increasing mass faster than the density of states in- 
creases. 



inverse mass density 




T=139.7 ± 3.2 MeV 
m„ = 500 MeV 



CHISQ = 1.19 



1.0 1.5 Z.O 

PARTICLE MASS (GeV/c»»2) 



FIG. 7: The solid curve is a fit to cross sections observed in 
this experiment using the inverse of the mass density p(m) 
suggested by the Hagedorn bootstrap model. 

The use of this expression for the cross section suggests 
that to some extent we may observe an exponential de- 
crease in da^ /dy with increasing mass because we choose 
to measure only one of the many available states at mass 
m, whereas the number of available states increases ex- 
ponentially with increasing mass. 



C. A Gas Model 

In this section we will explore the possibility that a 
relativistic gas model may describe the observed particle 
production cross sections. We make use of an expression 
for the number of particles per unit volume, dn, in the 
center of mass system given by [4l| 



dn I 



(18) 
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where T is the gas temperature and is a chemical po- 
tential. The volume element in momentum space, d^p, 
is given by 



d^p = 2-Kptdptdpi 
= 2t: ptEdptdy 



(19) 



The total energy of a particle is E. pt and pi are re- 
spectively components of its momentum transverse and 
parallel to the beam direction in the center of momentum 
system. The variable y — {l/2)ln[{E + pi)/{E — pi)] is 
the rapidity. In our experiment pi ~ 0. Thus one can use 
E = y/pf+rn? for the total energy of a particle of mass 
m and take y — Q. 

In order to calculate a primary cross section a simple 
assumption is made that 



daP 
dy 



A- 



y=o 



dn 
dy 



(20) 



To obtain the cross section for any particle it then be- 
comes necessary to evaluate the integral 



dn 
dy 



Pt E dpt 



(21) 



These integrals are done numerically. The second term 
in the denominator of the integrand should be -1 for 
bosons and -1-1 for fermions. In our applications we find 
that for masses above those of the kaon and pion the sec- 
ond term is insignificant, so that all of the calculations 
in this paper have actually used the boson form. 

A first attempt to describe the primary cross sections 
of Table [IV] was made by assuming the chemical potential 
^ = for all particles and by fitting equations [20] and [21] 
to constants A and T. The result is the curve shown in 

Fitted values are A = (3.69 ± 0.71) x 10^ mb/GeV'"^ 
and T = 0.1317±0.0038 GeV. Unfortunately the for 
this fit, = 14.15 for 4 degrees of freedom, is unaccept- 
able. The main difficulty seems to be a failure to predict 
a large enough pion cross section. This difficulty can be 
addressed by selecting a suitable chemical potential. 

From the thermodynamics of a system with k types of 
constituents, one knows that a change dU in the internal 
energy of the system can be written as (42(1 



dU = TdS - pdV + ^ n^dNi 



(22) 



It should be possible to estimate a value for the chem- 
ical potential of a pion, /^^, using a gedanken experiment 
which adds only one pion {dN^^ = 1) to a system with 
temperature T and pressure p while arranging for entropy 
S and volume V to remain constant (dS = dV = 0) . Then 
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FIG. 8: The curve is a fit of the statistical gas model in 
Ea. (|2ip to the primary cross sections per spin state of 6 par- 
ticles when the chemical potential ^ = for each particle. 



the chemical potential has a value = dU . This sug- 
gests that /iTT — fn-K , but the relation with the pion mass 
must be such that /i^r < Ti^r if the integrand of Eq. (|^T|) 
is to have no singular values. 

Provided the system already contains constituent u 
and d quarks, one can imagine that dU need not sup- 
ply the energy for existing quark masses. In that case we 
subtract existing quark masses [l^ as below. 



fJ-n 



0.1396- 0.010 
0.1296 GeV (23) 



We have used this same value of /i in Eg. lj^Tj) for each 
of the 6 measured primary cross sections in order to fit for 
new values of A and T. This fitted curve can be compared 
to the data in Fig. [5] 

The value of A = (1.054 ± 0.156) x 10^ mb/GeV^ and 
T = 0.1371 ± 0.0037 GeV. The = 4.35 for 4 degrees 
of freedom is quite satisfactory from a statistical point of 
view. 

The use of a common chemical potential /i for all the 
particles works best for fitting the pion data, and as we 
shall see later, it may influence the average kaon mo- 
mentum in a beneficial way. It is not possible with the 
present data set to exclude the possibility that fi = for 
the heavier particles. On the other hand, it is easy to rule 
out the use of the thermodynamic arguments like Eq. (|^^ 
and Eq. (j23p to determine a unique chemical potential for 
each mass because no reasonable fit to the cross sections 
was possible for that approach. 



X. AVERAGE TRANSVERSE MOMENTA 

Two significantly different temperatures are candidates 
for describing the production cross sections we have mea- 
sured, T=204 MeV used with a simple exponential in the 
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FIG. 9: The curve is a fit of the statistical gas model in 
Ea. (|21|l to primary cross sections per spin state when the 
chemical potential fj, = 129.6 MeV for each of the 6 particles. 



mass or T=137.1 MeV used in a relativistic gas model. 
We will first try to determine if either of these two tem- 
peratures is more appropriate for characterizing the mea- 
sured values of average transverse momentum. 

Consulting Eq. ((2T|) , we find that the relative transverse 
momentum spectrum at y=0 should be given by 



dn 



dydpt 



PtE 



y=o [e(-E-M)/r T 1] ■ 



(24) 



The average transverse momentum is then calculated 
for each mass value as 



<Pt> = 



/o°° Pt {dn/dy dpt ) dpt 
J^{dn/dydpt) dpt 



(25) 



An initial calculation was done using chemical poten- 
tials /i = 0. Curves of <pt> versus mass are plotted in 
Fig. [TO] for T==204 MeV and T=131.7 MeV, the latter 
being slightly more appropriate for /.t = as in the fit of 
Fig.m 

In Fig. [TO] it appears at first sight that the value of 
T=204 MeV has some amount of promise for being use- 
ful and that T=131.7 MeV is completely inappropriate. 
Neither curve recognizes the more pronounced decrease 
with mass in the experimental values of <pt> f^or the 
kaon and pion. 

This last feature is improved by using the chemical po- 
tential ^ = 129.6 MeV. Fig. [TT] compares use of T=204 
MeV and T= 137.1 MeV, the latter number being more 
relevant for the non zero value of /i as in Fig. [5] The 
higher temperature curve now comes closer to the mea- 
sured values for the kaon and pion, but it is still well 
outside experimental error bars. The lower temperature 
curve remains unacceptable because of a nearly uniform 
displacement below the data points. 
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FIG. 10: The curves represent average transverse momentum 
expected from a Bose-Einstein distribution with chemical po- 
tential = 0. The solid curve is for temperature T = 204 
MeV and the dashed curve is for T = 131.7 MeV. 
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FIG. 11: The curves represent average transverse momentum 
expected from a Bose-Einstein distribution with chemical po- 
tential n = 129.6 MeV. The solid curve is for temperature 
T = 204 MeV and the dashed curve is for T = 137.1 MeV. 



This displacement can be understood if one uses a sug- 
gestion made by Levai and Miiller [20] which said that 
part of the transverse momentum may be due to the 
transverse flow of the g whole. In order to cal- 

culate the effects of transverse flow it is advantageous to 
define a particle density 4- flow vector dN'^. Within a 
constant factor this is given by [4ll | 



p" 



dN^{p)=dMb^) 



1 



(26) 



where (p'^/p") is the ly velocity component of particles 
in the infinitesimal group diV". The particle distribution 
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function is here written in terms of an invariant p'^U^. 
In this form it is sometimes called the Jiittner distribu- 
tion [i^, Uu is the constant 4-flow velocity of the 
gas: Ui, = (jfc ,—jfVf), where Vf is the flow velocity 



and 



vl/c^. 



If the invariant is evaluated 



in the collision center of mass system for the case of a 
particle with momentum p parallel to the flow velocity 

Vf, 



p''U^ = jf{E-f3fpc)^E', 



(27) 



where E' is the particle energy in the rest system of the 
flowing gas. 

We obtain an expression which is needed for calculat- 
ing cross sections by forming the invariant product 



P^U, 



(28) 



When the flow velocity is zero, u/ = 0, the invariant 
product reduces to dn as used in Eq. ([2T|) . If the flow 
velocity is not zero and is parallel to the transverse mo- 
mentum pt, the invariant product becomes a more gen- 
eral expression for the invariant dn. 



dn = 



Pt E' dpt dy 



(29) 



E' is the energy of a particle in the rest system of the 
flowing gas and is calculated using Eq. (f27|) . 

The cross section data points are not effective in con- 
straining values of the flow velocity (3f = Vf/c. Because 
of this we have estimated a value of /?/ by fitting the 4 
best measurements of <pt>, those of the tt, K, p, and 
A. Nevertheless it is still necessary to select a reasonably 
correct value of T to be used in Eg. ips)) for computing 
<Pt>. We find a value of /3/ = 0.27 ± 0.02. The error is 
merely statistical and in no way reflects the systematic 
uncertainty caused by secondaries from resonance decays. 

Some insight into the systematic error might be gained 
by observing how close the fitted curve in Fig. [12] comes 
to the <pt> measured for the proton, since the heavy 
mass of the proton causes recoil effects in decays to be 
less likely to shift average values of momenta. The curve 
is well within the statistical error on the measured value 
of <pt> for the proton. 

The curve in Fig.[T2lwas obtained using T — 131.5±3.4 
MeV. This is the temperature obtained if one fits the pri- 
mary cross sections using Eg. ip^ with (3f — 0.27. Thus 
the cross section fit and the /3/ fit are consistent. For 
the cross section fit the normalizing constant fEq. ([20l) ) 
is A = 600 ± 86 mb/GeV^ and the = 3.57. The fitted 
curve using these values is imperceptibly different from 
that shown in Fig. [SI 
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FIG. 12: The solid curve represents average transverse mo- 
mentum expected from a Jiittner distribution with T — 131.5 
MeV and chemical potential ^ — 129.6 MeV when momenta 
are boosted by a transverse flow with l3f — 0.27. 



XI. TRANSVERSE MOMENTUM SPECTRA 

Since the experimental pt spectra of the particles are a 
mixture of primary particles and resonance decay prod- 
ucts, it is a challenge to find a useful means of comparing 
them to the statistical predictions for primary particle 
production. The curves in Fig. [T3] and Fig. [T3] are statis- 
tical model predictions for primary pt spectra of pions, 
kaons, and protons. Because these particle cross sections 
vary by an order of magnitude, the vertical scales for 
the curves have all been renormalized arbitrarily to fit in 
the same plot. Experimental data points have also been 
renormalized in order to help in making comparisons of 
data with model predictions. 

Figure [131 compares data points with a model that has 
T=204 MeV and a chemical potential = 129.6 MeV. 
Data points are cut off sharply at low momentum by ab- 
sorption in the spectrometer material. The kaon and pro- 
ton data are not extremely different from the predicted 
primary curves, but the pion data points are noticeably 
more condensed at low pt than this model predicts for 
the primary pions. 

Figure [HI compares data points with a model that has 
T=131.5 MeV, a chemical potential ^ = 129.6 MeV, and 
a transverse flow velocity /?/ — 0.27. All three particle 
spectra are similar to model predictions in an encourag- 
ing way. 

The influence of resonance decay products on the ex- 
perimental Pt spectra is perhaps more complex than one 
might initially guess, especially for the spectrum of the 
much lighter pions. We have examined spectra from the 
production and decay of two "typical" resonances in our 
spectrometer using a Monte Carlo program. The pro- 
gram generates particles at T=204 MeV using no chem- 
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FIG. 13: The curves represent transverse momentum distribu- 
tions for pions, kaons, and protons as expected from a Bose- 
Einstein distribution with temperature T = 204 MeV and 
chemical potential /i — 129.6 MeV. Normalization was arbi- 
trarily chosen to fit all three curves within the plot. Data 
points for the 3 particles were also normalized to facilitate 
comparisons with the 3 curves. 
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mass of M* = 1.6 GeV and a decay mode of M* p + ir. 
The histograms in Fig. [15] show the Monte Carlo spec- 
tra. The average pt values for the Monte Carlo proton 
and pion spectra are quite close to the experimentally 
measured values. 

The shape of the pion spectrum shows a surprising 
peak at low momentum. Much of this peak comes from 
M* resonances initially directed into the hemisphere op- 
posite our detector. Backward pions from the decay of 
this resonance pass through the aperture of the detector 
and are processed as very low momentum pions. A small 
bump seen at pt ~ 250 MeV/c in the experimental pion 
data might be associated with such an effect. 



M.C, PIONS 

FROM M*=l,3 GeV 




Pt (GeV/o) 



).5 1.0 
Pt (GeV/c) 

FIG. 15: Histograms are Monte Carlo generated distributions 
using a gas with T = 204 Mev. Dashed histogram represents 
pions from decay of M* = 1.3 GeV to vrp. Solid histogram 
represents protons from decay of M* = 1.6 GeV to np. Over- 
all normalization is arbitrary. Data points are normalized to 
give areas comparable to the histograms for comparison pur- 
poses. 



XII. LUMINOSITY AND TRIGGER ISSUES 



FIG. 14: The curves represent transverse momentum distri- 
butions for pions, kaons, and protons as expected from a 
Bose-Einstein distribution with temperature T — 131.5 MeV, 
chemical potential fj, — 129.6 MeV, and transverse flow ve- 
locity /3/ = 0.27. Normalization was arbitrarily chosen to fit 
all three curves within the plot. Data points for the 3 parti- 
cles were also normalized to facilitate comparisons with the 3 
curves. 



ical potential and assumes a flat rapidity distribution. 
Experimental acceptance cuts and magnetic field are in- 
cluded in the analysis, but absorption effects are ignored. 
For pion sources the typical resonance is chosen to have a 
mass M* = 1.3 GeV and a decay mode M* ^ p + n. For 
proton sources the typical resonance is assumed to have a 



There are two methods used to calculate the luminos- 
ity of a bunch crossing. One is the "machine method" 
and the other is the "cross section method" . The ma- 
chine method consists of first measuring and estimating 
properties of the particle bunches and the particle accel- 
erator and then using these to calculate luminosity from 
first principles [11]. The cross section method counts the 
number of events from a beam crossing that result from 
a particular process and uses the known cross section for 
that process in Eq.([T]) in order to calculate the luminos- 
ity. Once the luminosity is determined, all other cross 
section measurements are reduced to event counting ex- 
periments. 

When we first began to evaluate cross sections for the 
1988-1989 Tevatron run, the most accessible luminosities 
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were machine calculated values for which the data were 
posted at approximately two hour intervals by the ac- 
celerator operators. We have used these same machine 
luminosities [IT] to calculate all cross sections presented 
in this paper. 

The trigger used to collect the data which we present 
from this experiment was designed to select non- 
diffractive (nd) events. Therefore a cross section com- 
puted for this trigger, atrig, should closely approximate 
the best known value for the non-diffractive cross section, 

Our primary trigger (PT) was typical of those used 
with collider detectors. It required a triple coincidence 
between far upstream and far downstream counters (BB) 
and the bunch crossing time TO. This coincidence was 
supplemented by time based vetoes to eliminate beam- 
gas events from halo and satellite bunches. Events se- 
lected by the PT trigger are called non-single-diffractive 
(nsd), because elastic (el) and single-diffractive (sd) 
events are effectively elimnated by it. Since our highest 
luminosity was around 2 x 10^® cm~^sec~^ (about two 
orders of magnitude below that of CDF), an interaction 
occured in about every 200 bunch crossings, and there 
was little need to consider cases where two sd events in 
the same bunch crossing could fake one nsd event. 

The trigger (SI) used for collecting data reported in 
this experiment made use of the PT trigger described 
above in coincidence with a "track" in the spectrometer 
arm. A track was defined by requiring drift chamber hits 
in at least 3 of 4 possible planes both before and after 
the analysis magnet. This requirement of a spectrome- 
ter track near zero rapidity effectively eliminated double- 
diffractive (dd) events from the event trigger (SI). 



TABLE VII: This table contains cross sections used to check 
our machine-derived luminosity against the more commonly 
used cross-section-derived luminosities. Column 3 gives 
sources for the cross sections listed in column 2. 



Cross Section 


Value (mb) 


Source 


Cytrig 


41 ±6 


m 


Cdd 


4.43 ± 1.18 


[45J 




9.46 ± 0.44 


[46J 


^ ineL 


59.3 ±2.3 


[47] 


O'nd 


45.4 ±2.6 


[45],[46],[47] 



The first row of Table IVlTl gives the cross section of the 
SI trigger as computed using machine-determined lumi- 
nosity for the CO intersection, atrig = 41 ± 6 mb. We 
should compare atrig to the non-diffractive cross section 
and a-s derived from cross-section-determined luminosi- 
ties given in Table IVIII 



cTdd - cTsd = 45.4 ± 2.6 mb. 



(30) 



In this calculation we have used an average value for ainei 
which was adopted by the CDF and DO collaborations for 
defining cross-section-determined luminosities [47| . Our 
machine-determined atrig is about 10% lower than the 



cross-section-determined and, but the two values agree 
within the errors. Should one ever choose to compare 
cross sections in this paper with those published by CDF 
or DO, then it might be desirable to shift one set or the 
other systematically by 10%. 

It may be interesting to note that increasing all of our 
absolute cross sections by 10% would make the scale fac- 
tor for our cross section fit agree almost exactly with that 
for TTp scattering in Eq.([7]). 



XIII. ENERGY DEPENDENCE 

Much can be learned about the properties of the ther- 
mal state by considering the production properties of the 
(j) meson. As one can see from Fig. [3ljd), there are neg- 
ligible contributions to the inclusive 4> cross section by 
secondaries from higher mass decays. The (j> cross sec- 
tion at y=0 is essentially a primary cross section. If the 
same thermal production mechanism prevails in other ex- 
periments, then a measurement of the phi cross section 
in those experiments will allow a comparison of scales 
among experiments. The decay of the (j) into lepton 
pairs and kaon pairs means that most experiments are 
equipped with a simple trigger which can record its pres- 
ence. 

Absolute cross sections at collider energies are rarely 
reported, but we are able to plot two additional values of 
da/dy\y^Q for (j) mesons in Fig. [H] as a function of s 48] . 
The lower energy cross sections are for pp collisions, and 
the highest energy one is our value at -y/s = 1800 GeV 
for pp collisions. 

We have observed evidence for the thermal state only 
in collisions with the potential for annihilation {pp and 
e~^e~ initial states) . The similarity of the pp and pp cross 
sections for mesons shown in Fig. [11] suggests that the 
origin of the thermal state may be in the interaction of 
gluons, which are common components of both protons 
and antiprotons. 

An inelastic cross section such as the total cross section 
for (j) production might be expected to have an energy 
dependence which is bounded by some relation like the 
Froissart limit [i^, atot{4') In^ s. To the extent that 
inclusive cross sections have a flat rapidity distribution 
for \y\ < Umax, we can write 



, , da 

<7tot{9) — 2 Umax-r 

ay 



y=o 



w / 2\da 
ln(s/m ) — 
dy 



y=o 



In this case it seems likely the </) cross section at y=0 
would respect a reduced Froissart limit of da/dy\y^o < 
In s. 

In Fig. [16] we have plotted best fits to 



da 
dy 



constant 



0.58 ±0.07 mb) 



and 

da 

dy 



y=o 



A ln(s/so) (= 0.039±0.005)-(ln(s/TO^) mb), 
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FIG. 16: Cross section da/dy at y=0 for inclusive meson 
production is plotted vs. s, the square of collision energy. 
Cross sections show very little dependence on either energy 
or an annihilation option. 

where sq is an energy scale factor to which fits are insen- 
sitive. We have chosen sq to be the square of the pion 
mass. 

The Ins fit has a, = 0.59, and the constant cross 
section fit has a — 1.1. Both possibilities are accept- 
able for 2 degrees of freedom. The data are consistent 
with a constant (f> cross section or a slight increase with 
energy near a Froissart bound. 

XIV. NUCLEAR CROSS SECTIONS 

Although one can see in Fig. [4] that the two nuclear 
cross sections do not fit the thermal model well, they 
seem to exhibit a tendency to be in the neighborhood 
of the dashed curve that describes the model. The cross 
section per spin state for deuterium is a factor of 5.0 be- 
low the dashed curve, and the cross section per spin state 
for tritium is a factor of 8.7 above the dashed curve. It 
should be mentioned that these plotted points are esti- 
mates of the cross sections and not direct measurements. 
We should ask if the estimate for the deuterium cross 
section could be larger or if there might be a reason for 
the tritium cross section to appear so large. 

The invariant cross sections, {l/n)d'^a/dydp^, for d+, 
d~ , and t'^ were measured directly in a very limited 
range of pt, but values for the integrated cross section 
da/dyly^Q had to be obtained without benefit of fully 
measured distributions of transverse momentum. The 
published estimates for deuterium and tritium were made 
using variations of a Boltzman momentum distribution 
which had adequately fit pt spectra for K~ and p pro- 
duced at y/s = 1.8 TeV. Systematic errors were assigned 
based on the results of these variations. 

The tritium cross section was based on 8 events sat- 



isfying experimental cuts. Only one possible anti tri- 
tium event was observed. Anti tritium is more heavily 
absorbed by the spectrometer material, but the record- 
ing of events was heavily prescaled by total multiplicity, 
sometimes as much as a factor of 26, so at this level of 
statistics it is only suspicious that there are not more anti 
tritium event candidates observed. 

There are two reasons to consider why the tritium cross 
section plotted in Fig. [4] might be higher than expected. 
The first reason is that the observed tritium nuclei might 
not originate from the pp collisions. They might come 
instead from secondary interactions in the wall of the 2 
mm thick aluminum beam pipe. 

Approximately 2 x 10^ particles entered the spectrome- 
ter arm during the experiment. Some of these in principle 
could have ejected a tritium nucleus from aluminum in a 
quasi elastic process such as 

n+ +'^'^Al^TT+ +^H+^^Mg, (31) 

where '^'^Mg is a spectator nuclear fragment. Using only 
this process, we find ^ 12% of the spectrometer particles 
have momentum high enough to produce the '^H candi- 
dates we observe. 

Spallation experiments have reported total cross sec- 
tions of 27 mb for producing ^^Mg using protons with 
momentum of 1.2 GeV/c (50]. In no case have we found 
a published rate for '^H production in a spallation ex- 
periment, although some might have had the ability to 
observe it. Much of the production of light masses is 
expected to come from a low energy evaporation process 
after the initial impact. However, if all of the 27 mb cross 
section for producing '^'^Mg proceeds by Eg. ([51]) with 
in the forward direction, then there would be 800 tritium 
candidates in the spectrometer. If fragments exist in 
such numbers, it cannot be safely assumed that all of 
these would be eliminated by the vertex cuts we used. 

There is a second reason the tritium point in Fig. 3] 
might appear too high. It could be that there are more 
contributing states to the cross section than the (2J+ 1) 
spin states we assumed. The estimated cross section was 
divided by only (2 J -I- 1) spin states before plotting. 

Some evidence exists that well defined "nuclear" states 
exist in which one of the usual nucleons is replaced by a 
particle resonance. The best established examples of this 
can be found in the A° hyperfragments, where the A° 
lives long enough to give a distinctive decay signature in 
a detector. A collection of old data in reference [5l| fists 
12 hypernuclei ranging from ^Hp^ to ^^C\ and tabulates 
their binding energies along with some observed decay 
modes. More recentlly high resolution spectrometers |52l | 
have resolved excited energy levels for the bound A° in 
heavier nuclei. For ^^1a energy levels for s, p, d, and f 
A° orbitals have been resolved. 

If N* resonances can also be bound in nuclei, there 
might be many modes which would contribute to the ob- 
servation of '^H production in pp collisions. Although 
there seem to be no experiments designed to explore this 
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notion, one can find some indication of it. Inelastic cross 
sections for pions on tritium peak at a pion energy around 
140 MeV [53], suggesting the formation of A(1232), and 
final state particle ratios in tt^ collisions with and 
■^iJe are consistent with dominance by the A baryon res- 
onance 1541. 



a larger radius r > r^. According to Eq. (j33p a larger ra- 
dius will imply a correspondingly lower temperature and 
a smaller cross section than does the smaller radius of a 
hadron having the same mass as the nucleus. 

Assuming a nuclear volume increases as the number of 
constituent hadrons, A, we can write an effective nuclear 
radius ta as 



XV. LITTLE BANG NUCLEOSYNTHESIS 

The Little Bang collision as observed in this experi- 
ment takes place in approximately 10^^^ seconds as the 
Lorentz contracted disks of the proton and antiproton 
pass through each other. Interfering pions continue to 
be emitted [2l[ from the vicinity of the collision for 
times ^ 10'^ longer than the collision time, since the 
measured size of the emission volume has a "radius" of 
approximately 1 fm, somewhat larger than this in the 
longitudinal direction than in the transverse direction. 
Slightly larger emission dimensions are correlated with 
larger multiplicities of produced hadrons. 

Several estimates [2l|, [sll have been made of the energy 
density in the emission volume at the time of hadroniza- 
tion. The values reported are sensitive both to the mod- 
els used and the fraction of observed final state particles 
assumed to originate from the emission volume. Energy 
densities obtained range from 4 to ~ 60 times that of 
a proton. 

Primary cross sections per spin state for the light 
hadrons decrease with mass according to a Boltzmann 
distribution having a characteristic temperture T = 204 
MeV. See Fig. 21 We can understand this temperature in 
a qualitative way if we assume that the collision deposits 
energy into a volume of hadronic size at the limit of the 
Uncertainty Principle. In that case Ai? cAi = Tic, where 
/S.E is the energy deposited and cAi is a characteristic 
hadronic radius ru- 



AE = hc/vh 



(32) 



In the Uncertainty Principle AE is the measure of the 
drop-off energy of a probability function. For a Boltz- 
mann distribution we can identify AE with T [s^ . 



T = hc/rh 



(33) 



If all produced hadrons share the same radius r^, then 
T is a constant which can be found experimentally. A 
value of T = 204 ± 14 MeV implies a hadronic radius 
rh = 0.97 ± 0.07 fm. In some sense has to be an 
average radius. 

In this experiment we observe the production of deu- 
terium and anti deuterium nuclei. The radius of a nu- 
cleus does not fit comfortably inside a hadronic radius, 
so its production must wait for the volume to expand to 



(34) 



The temperature for nuclear production becomes Ta — 
hc/rA, and the cross section in millibarns per spin state 
for producing a nucleus A becomes 



da (A) 



dy 



y=o 



45.3 e" 



ia/Ta 



(35) 



Normalization at 45.3 mb fixes the primary proton cross 
section at the value previously used with T = 204 MeV. 

A dashed curve in Fig. [T7| plots the locus of nuclear 
cross sections per contributing state as predicted by 
Eq. ([55)1 . The experimental estimates for the d and d cross 
sections agree fairly well with the prediction of a modified 
temperature for nuclear production, whereas the experi- 
mental estimate of the ^iJ cross section becomes 4 x 10'^ 
times larger than the dashed curve. 
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FIG. 17: Primary cross sections per spin state are compared 
to a Boltzmann distribution in mass for a temperature of 
T — 204 Mev (solid curve) and to a distribution modified 
to accomodate the radius of each individual nucleus (dashed 
curve) . 

Nuclear cross sections decrease so rapidly with mass 
according to Eg. pS)) that only the lightest nuclei might 
be detected at colliders. Table I Villi lists predicted pro- 
duction rates relative to p rates for several light ele- 
ments. Based on our observation of 19 anti deuteron 
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events at J Cdt = 5.52 ± 0.81 nb we estimate that at 
/ £dt = 2100p6"i the CDF detector at the Tevatron 
could observe ~ 3.6 x 10^ anti deuteron events and ~ 27 
anti '^He events, provided there is no prescaUng of the 
event triggers such as was used in our experiment. This 
estimate ahows for the increased geometrical acceptance 
of the CDF detector but not for other detector specific 
features. 



TABLE VIII: The production rate of a light nucleus relative 
to the production rate of antiprotons is calculated using the 
temperature model of Eq. (|35[) . J a is the nuclear spin, and 
A is the mass number of the nucleus. Spin degeneracies are 
included in the calculation of the rate ratio R. 



Nucleus A Ja Na /Np 

Hi 1 172 To 

2 1 1.38 X 10"^ 

^H, ^He 3 1/2 2.27 x 10"'' 

^He 4 1.04 X 10"" 

^Li 6 1 2.52 X 10"^° 

''Li 7 3/2 3.62 x 10"^^ 




2 4 6 



PARTICLE MASS (GeV/c^) 

FIG. 18: Plotted points are relic cosmic abundances of ^H, 
"^He, *ffe, and ^ Li relative to hydrogen 57] after dividing 
each ratio by (2 J + 1). Variations of the Big Bang Nucle- 
osynthesis model strive to reproduce these data points. The 
solid curve when used with AM = niA ~ rup has the slope 
for T — 204 MeV, which fits elementary particle production 
in this paper. The dashed curve represents a temperature 
modification which might accomodate nuclear size. 



If nuclei are reaUy formed in nucleon-nucleon coUisions, 
it should be relatively easy to verify the fact with suitable 
triggers in conjunction with one of the high resolution 
detectors that currently exist. It would be essential to 
require that the instrumentation not saturate on pulse 
heights up to 10X-20X minimum ionization. The high 
statistics expected for the lightest nuclei (^H, , ^He) 
would permit one to study the extent to which they are 
produced in pairs or to put limits on the equality of nu- 
clear and antinuclear cross sections. Because a tempera- 
ture of T = 204 Mev also seems to organize the hadron 
production cross sections from e^e~ collisions, it may be 
useful to search for nuclear production in those data sets 
as well. 

Little Bang Nucleosynthesis (LBN) would seem to be 
quite different from Big Bang Nuleosynthesis (BBN) [H^] . 
BBN begins with a large net baryon number, but LBN 
begins with zero net baryon number, initiated perhaps by 
annihilation or by gluon-gluon interaction. Electrons and 
photons participate in BBN in a large way, but there is 
little evidence for either of these in LBN (see section IV) . 
In BBN nuclei are formed at temperatures comparable to 
nuclear binding energies. Nuclei from LBN appear to be 
formed somehow in an environment with temperatures 
that are many times larger than nuclear binding energies. 

Fig. HH] has been included here to illustrate the extent 
to which LBN fails to reproduce the relic cosmic abun- 
dances that BBN accounts for so well [13] ■ The ratio 
of cosmic abundances relative to hydrogen have been di- 
vided by (2 J-t-1) for plotting. Nuclear curves from Fig.fTTl 
have been inserted into the figure for comparison. 



XVI. CONCLUSIONS 

This paper has presented 9 absolute cross sections 
(da/dyly^o) for particle production at 90^^ to colliding 
p and p beams with ^/s = 1.8 TeV. One of these cross 
sections is for photon production. Two are for nuclear 
production. The remaining six arc for light elementary 
hadron production. 

The six hadron cross sections were corrected for sec- 
ondary contributions from decays of higher masses in or- 
der to produce a set of six primary cross sections. When 
the primary cross sections are divided by the spin de- 
generacy (2 J -|- 1) of individual particles, we learn that 
the primary cross sections per spin state precisely follow 
an exponentially decreasing distribution in particle rest 
mass with a temperature T — 204 ± 14 McV. By sum- 
ming over all spin states, one finds the total primary cross 
section for a particle of specified charge, spin J, and rest 
mass m is given by 



dcrP(tot) 



dy 



= 0.721 •(7rA^)(2J + l)- e" 



(36) 



Using an argument based on the Uncertainty Principle, 
we set T = Tic/rh^ where rt is a hadronic radius given by 
r/i = 0.97 ± 0.07 fm. For nuclear production we replace 
fh by rA — I'hA^^^ and find a temperature Ta- When 
Ta is used in Eq. ([36l) . it yields an acceptable value for 
the anti deuteron cross section. 

This purely thermal production mechanism exhibits no 
significant strangeness suppression, implying that exper- 
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imentally observed strangeness suppression is entirely ex- 
plained by contributions from decay products of higher 
masses. As discussed in section IV, we were unable to ob- 
serve any primary production at the photon and electron 
masses, so the thermal mechanism is probably relevant 
only for strongly interacting particles. 

It appears there may be two temperatures that are 
useful in describing the entire particle production pro- 
cess. The temperature T=204 MeV precisely describes 
the primary particle ratios in the initially produced state, 
and a temperature T=132 MeV describes the momentum 
distributions in a final decay state in which pious have 
a uniquely low momentum. An exothermic transition 
between the production state and decay state might ac- 
count for a transverse flow velocity oi (3/ — 0.27. The 
average kinetic energy per particle associated with this 
transverse flow is ~ 20 MeV, and thus it is not a trivial 
amount. 

It is possible to observe the thermal production mecha- 
nism in pp collisions at 90*^ (y=0) because of its dominant 
cross sections in the low mass region. For the production 
of masses in or above the charm sector there is evidence 
that other mechanisms such as parton scattering become 
dominant. The dominance may extend to many orders 
of magnitude at the highest particle masses. 

We have examined cross sections for hadron production 
in the charm sector by e^e~ collisions using published 
values of multiplicity fractions. To the degree that decay 
branching ratios are well known, we find that primary 
e~^e~ multiplicity fractions also follow a Boltzmann dis- 
tribution in mass for a temperature T — 204 MeV. There 
is no significant evidence for strangeness suppression in 
the charm sector examined. 

Since the 4> meson is essentially a primary particle with 
negligible secondary contributions, one can possibly gain 
some insight into the origin of the thermal mechanism 
by studying </> production. The ip cross section at 90*^ is 
consistent with being the same for pp and pp collisions 
over a range of center of mass energies that vary by a 
factor of ~ 30. This may suggest that gluon interactions 
rather than annihilations are the essential source of the 
thermal production. 
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XVIII. APPENDIX 



Cross Section Estimates 



We have used partial cross sections and raw data found 
in the literature to make estimates of total inclusive cross 
sections for four massive particles. The fact that exper- 
imenters themselves did not present total cross sections 
testifies to the difficulty and large uncertainty associated 
with such calculations. Nevertheless the estimates made 
below may be useful within an order of magnitude. 

J/ip'. Using a minimum bias trigger in a special run, 
CDF measured a cross section [2J| for the prompt pro- 
duction of J/tp with a value of 

a{pt > 1.25, \y\ < 0.6) = 2.86 ± OMtait (37) 

Prompt J/ijj's are those which are consistent with orig- 
inating at the beam line rather than from detached ver- 
tices. Since the entire pt spectrum was measured for the 
J/ip in this experiment, we used that information to cor- 
rect the measured value upward by a factor of 1.28 before 
dividing by the effective interval Ay = 1.2. The result 
given below was divided by (2 J -I- 1) = 3 and plotted in 
Fig.El 



da{J/ilj) 



dy 



y=o 



3.06 ^b. 



(38) 



D*+: Using a detached vertex trigger, CDF measured 
the D*~^ cross section as a function of pt above pt = 6.0 
GeV/c 24]. We have used the J/ip momentum spec- 
trum as a guide to increase this cross section by a factor 
of 30.6 to account for the unmeasured portion of the pt 
spectrum. Dividing by the effective rapidity interval of 
Ay = 2, we estimate 



da{D*+) 



dy 



0.075 mb. 



(39) 



This value is divided by (2 J -I- 1) = 3 and plotted in 
Fig.El 

: CDF has reported 17.5 of these events with a 
decay mode of Sjj" J/tpE~ [ssj . 

We divide this number by integrated luminosity L ~ 
1.9 fb~^ and multiply by 10 to correct for J/ijj acceptance 
over the rapidity range when J/^ momentum is pt > 3.0 
GeV/c. We multiply by 3.26 to correct the unobserved 
J/ip momentum spectrum below pt = 3.0 GeV/c using 
the Pt spectra of reference |23| as a guide. We divide 
by 0.0593 to correct for the branching fraction of J/ip 
into the observed pair [18] and multiply by 10^ to 

estimate a typical probability for hadron decay into J/ip. 
Dividing the resulting cross section by Ay = 2, we obtain 



da(S,-) 



dy 



2.53 X 10"^ mb. 



We have divided this number by (2 J + 1) 
plotted it in Fig. \5\ 



(40) 
2 and 
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CDF has measured the cross section at -^s = 1.96 
TeV for pp — > H^X , where X is any thing and Hi, includes 
b hadrons and anti b hadrons 12311. 



a{pp HbX, \y\ < 0.6) = 17.6 ± 0.4 



^li Mb. 



(41) 



We divide this by 2 to obtain a cross section only for 
hadrons and by Ay = 1.2 to estimate da/dy. The J5+ 
cross section can be estimated using the branching frac- 



tion B{b B+) = 0.398±0.012 [ll| as a multiplier. The 
resulting cross section estimate for production is 



da{B+) 



dy 



2.92 Aib. 



(42) 



We have further reduced this value by 10% in an effort 
to scale it to a collision energy of y^s = 1.8 TeV and 
divided by (2 J + 1) = 1 for plotting in Fig. [5l 
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